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; Abstract An analysis of all available data (CELLO, CLEO, BABM) in the range [1 40] GeV^ for 

the pion-photon transition form factor in terms of light-cone sum rules with next-to-leading-order ac- 
00 ' curacy is discussed, including twist-four contributions and next-to-next-to-leading order and twist-six 

, corrections — the latter two via uncertainties. The antithetic trend between the BABAR data for the 

7*77r" and those for the 7*7?7(?7') transition is pointed out, emphasizing the underlying antagonis- 
' (— I tic mechanisms: endpoint enhancement for the first and endpoint-suppression for the second — each 
^ i' associated with pseudoscalar meson distribution amplitudes with distinct endpoint characteristics. 

i' Keywords Transition form factors • pion distribution amplitude • light-cone sum rules 

1 Introduction 

m 

\ This is a report about the pion-photon transition form factor — an observable that is considered to 

be the prime example of an exclusive process in QCD because in leading order (LO) the partonic 
interactions are purely electrodynamic, with the strong interactions being factorized out into the pion 
. distribution amplitude (DA). While the accuracy of the CELLO data [1| was not sufficient to provide a 

proof-of-concept of the QCD framework, the appearance of the CLEO data 01 established the validity 
of the collinear factorization of this process [3]. This situation changed dramatically in the year 2009 
when the BABAR Collaboration published new data [4] with very small error bars at intermediate 
momentum transfers reaching values just below 40 GeV^ — unreachable in the past — with still tolerable 
errors. These data for the scaled form factor Q^F^ '^'^ (Q^) exhibit two distinctive features: (i) The 
earlier data taken by CLEO at the same momentum-transfer values were confirmed, but with a higher 
accuracy, (ii) Unexpectedly, the data starting above 9 GeV^ and extending up to about 40 GeV^ do not 
$H ' scale with . Instead, they show a clear general tendency to rapidly grow with Q^, with the exception 

5^ , of two single data points at about 14 GeV^ and 27 GeV^ that lie just below the asymptotic prediction 

of perturbative QCD Q: limQ2^o^ Q^F"*" — !• ^/2f.^. Understanding the underlying enhancement 
mechanism(s), responsible for this behavior, is of significant importance because an increasing form 
factor challenges the basic concepts of QCD, like collinear factorization. The recent analysis in Ref. 
[sf, reported here, builds on light-cone sum rules (LCSR)s by taking into account various contributions 
from perturbative and nonperturbative QCD and by utilizing in the analysis data sets from different 
experiments [l|; 0; 0] for momentum-transfer values ranging from 1 to 40 GeV'^. Emphasis is given 
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to the statistical features and the way they depend on the various theoretical parameters. The main 
goal is the confrontation of theoretical predictions, obtained with the help of light-cone sum rules 
(LCSR)s, with all available experimental data for the pion-photon transition form factor. In addition, 
the results are compared with the BABAR data [1] for the rj- and ry'-photon transition form factors 
making use of the description of the rj — r/' mixing in the quark flavor basis 0] that involves the state 
\n) = {\uu) + \dd))/^/2. This makes it possible to link the 7*7 — >■ |n) transition form factor, multiplied 
by 3/5, to the form factor 7*7 — > 7r°, where the prefactor arises from the quark charges. 



2 Computational Method 

The method of light-cone sum rules, developed in 8^ for the pion-photon transition form factor, but 
not only , has been refined by years of experience in many analyses [13; [HI; [H; [H; [1| . It contains 
a set of rules that encode a computational model and as any model it encodes implicit predictions 
that turn explicit under some justified assumptions about its inherent (and auxiliary) parameters. The 
subtlety is how to adjust these parameters in such a way as to achieve the best possible accuracy of 
the obtained predictions with respect to the experimental data while avoiding sensitivity on particular 
values. 

The modus operandi of the LCSRs for the pion-photon transition form factor is expressed by [3 



Q2 /m|-_Q2^\ n ^ dx 



(1) 



where the abbreviations s = xQ'^/x and xg = Q'^/{Q'^ + so) have been used, and with the spectral 
density p{Q'^,x) = {Q^ + s) p^^ {Q'^ , s) taken with next-to-leading-order (NLO) accuracy. The hadronic 
threshold in the vector-meson channel has the value sq = 1.5 GeV^, M is the Borel parameter, and 
nip =0.77 GeV denotes the physical mass of the p meson. The first term in ([T]) stems from the hadronic 
content of a quasireal photon at low s < sq, while the second one resembles its pointlike behavior at 
higher s > sq. The main advantage of employing LCSRs is that one starts with the situation where 
both photon virtualities are sufficiently large, so that perturbation theory is safely applicable, and 
approaches the asymmetric kinematics with fixed and large and via a dispersion relation: 

F-*-*-(g^q^)= Hds^i^, (2) 

where the physical spectral density p{Q^,s) approaches at large s the perturbative one: p^'^{Q^,s) = 
^ImF'^ ''"^ (Q^, —s — ie) . The form factor {Q^,q^) can be calculated via perturbative QCD 

having recourse to collinear factorization, meaning that it can be written in the convolution form [l^ 

F^'''^^Q',q') = J\xTiQ^,q',^il,x)^i'\x,^il) + O {S'/Q') , (3) 

where the pion decay constant is /^r — 132 MeV and is the twist-four coupling. All nonperturbative 
information of the bound state is encapsulated in the coefficients a„ that are incalculable within 
perturbation theory and have to be modeled, or be computed on the lattice, e.g., [Tsl : [T6j |. via the 
moments {^^)-k = dx{2x — \)^ if'-n\x , p?') with the normalization condition dxip^\x, p^) = 1. 
The leading twist-two pion DA fulfills an evolution equation [iJl and can be expressed in terms of 
the Gegenbauer polynomials Cn^^{2x — 1) to read 



^^^Hx,p')^^^ix)+Y^ a,,{p')M^), (4) 



where '>pnix) = 6x{l — x)Cn^^{2x — 1) and ip'^^^x) = 'ipoix) = 6x{l — x) is the asymptotic pion DA [sl:!!^. 
The cornerstones of the analysis in [^, on which this report is based, are: 

— NLO radiative corrections: They are taken into account in the spectral density [l^, after correcting 
the error pointed out in [13*1, and using p^^^^ = p^^^ — Q'^ to avoid (large) logarithms of these 
scales. They are negative leading to suppression of the transition form factor. 
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Various contributions 



Fig. 1 (color online) Illustration of the various contributions to the pion-photon transition form factor. 



— Twist-four contributions: They are included via the coupling £ [0.15 0.23] GeV^ and per- 
forming one-loop evolution assuming the asymptotic form for the twist-four pion DA: tp^\x, /x^) = 
(80/3)(5^(/i^) x^{l~x)'^ . [Nonasymptotic forms were found before to be of minor importance.] These 
contributions are also negative and yield suppression. 

— Evolution of the coefficients a„: They are included at the NLO level of accuracy, using the QCD 
scale parameters A^^q^^'' — 370 MeV and ^qctT*' = 304 MeV, and provide suppression. 

— NNL0/3(, (for next-to-next-to- leading) radiative corrections [l3] and twist-six contributions [Tsl ]: 
They are taken into account implicitly in terms of theoretical uncertainties. This is justified because 
for the average value of M^(Q^) ~ 0.75 GeV^, used in our analysis Q, the net result of adding these 
(not large) terms is small, decreasing with from -f0.005 at = 1 GeV^ — where the positive 
twist-six term prevails — down to —0.003 at = 40 GeV^ — where the negative NNLO correction 
becomes stronger. Note that for AP — 1.5 GeV^ [l^, the twist-six term becomes negligible, while 
NNLO^Q still provides suppression. 

— Inclusion of resonances: A finite- width Breit-Wigner form is adopted in the spectral density in order 
to resolve the p and w resonances [1; [l3] • This entails a small enhancement of the transition form 
factor as compared to the ^-function ansatz. 

Figure [T] illustrates the approximate composition of the pion-photon transition form factor. 



3 Theoretical Predictions and Comparison with Data 

We now have all the elements for the data analysis. In order to get the most output from the data for 
the pion-photon transition form factor [1; 2; 4] and achieve the best insight about the particular role 
and relative strength of the first three Gegenbauer coefficients 02,04,05, we divide the data into two 
sets: (a) [1 9] GeV^— termed 'GLEO regime' and (b) [1 40] GeV^— 'entire range'. The results of 
the statistical 3-D analysis of these data sets are displayed in Fig. [5]— upper part — in the left and right 
panel, respectively, using the central value = 0.19 GeV^ of the twist-four coupling. Both panels 
show the corresponding la ellipsoids which are characterized by their principal axes and ellipses. The 
projection of the la ellipsoid on the plane (02,04) is represented, in both panels, by the larger ellipse 
enclosed by a solid line in red color. The smaller enclosed ellipse in the left panel denotes the cross 
section of the ellipsoid with the (02, 04) plane, whereas the shaded (green) rectangle shows the region 
in the (02,04) plane allowed by nonlocal QGD sum rules [isj . with the point in the center (marked 
by K) denoting the Bakulev-Mikhailov-Stefanis (BMS) pion DA The profile of this DA exhibits 
a double- humped structure with suppressed endpoints a; = 0, 1 due to the use of nonlocal condensates 
[Tsl — see also [l^ for a discussion of the endpoint region of the pion DA. All results are shown at the 
scale /igy = (2.4 GeV)^ after NLO evolution. Performing the analysis with the second data set (upper 
right panel), entails a sizeable coefRcient og, so that the la ellipsoid is lifted off the plane (02,04), 
while its projection remains almost the same as in the left panel. Remarkably, it still overlaps with the 
shaded rectangle containing the BMS model and the constraints from nonlocal QGD sum rules ^18| . 
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Fig. 2 (color online). Upper panels show 3-D graphics of la error ellipsoids in Gegenbauer space (02,04,06)- 
Lower panels show la error ellipses in the (02,04) space resulting from unifying la ellipses pertaining to the 
values of the twist-four coupling 5^ = 0.15, 0.19, 0.23 GeV^. Left panels refer to the analysis of all data 
in the range [1 -^ 9] GeV^, whereas the right panels give the analogous results for the region [1 40] GeV'^. 



Attempting to understand more quantitatively how the twist-four contribution affects our results, 
we concentrate on the (02,04) plane and present our findings in the lower part of Fig. [5] in terms of 
Icr ellipses. These are calculated by varying the 5'^ coupling from 0.15 to 0.19 to 0.23 GeV^ and then 
merging the obtained la error ellipses together to form the distorted ellipses shown in this graphics. Left 
panel ('CLEO regime'): The largest ellipse — dashed (red) line — results from combining the projections 
on the plane (02,04) of the 3-D data analysis discussed before. The smaller ellipse (solid blue line) 
shows the outcome of a 2-D analysis by means of 02 and 04. Its middle point with the coordinates 
(0.15,-0.09) and Xndf ~ 0-^ almost coincides with the center of the rectangle from [18]. Finally, the 
smallest ellipse (thick solid line), entirely enclosed by the previous one, is obtained by unifying the 
intersections with the (02, 04) plane of all ellipsoids generated by the variation around the central value 
of = 0.19 GeV'^. Right panel ('entire range'): Including into the statistical analysis the high-Q^ tail 
of the BABAR data modifies this picture completely. Now the composed error ellipse resulting from the 
2-D analysis (solid blue line) moves out of the region of the negative values of 04 — characteristic of 
the BMS models [3 — entering the positive domain. This entails a significantly worse value Xndf ~ 
relative to the value x^^f ~ 0.5 for the 'CLEO-regime' data set. On the other hand, the unified la error 
ellipse of the 3-D projections on the (02,04) plane (larger dashed red ellipse) keeps its position fixed still 
enclosing the major part of the 02 , 04 values obtained from nonlocal QCD sum rules (shaded rectangle 
in green color). Confronting these results with lattice estimates reveals that the 3-D projected error 
ellipse lies almost entirely within the boundaries from jisl (dashed vertical lines), even intersecting for 
the larger values of 02 with the narrower interval determined in [l^ (dashed-dotted vertical lines). In 
contrast to the left panel, the ellipse of the 2-D analysis for the entire set of data only poorly complies 
with the small 02 window of [l6j ] , while partly overlapping with the low end of the 02 range determined 
in [Ts"] . For comparison, some characteristic pion DAs are also shown in Fig. [5J asymptotic DA (♦) ; 
Chernyak-Zhitnitsky (CZ) model (■); projection of Model III from jJi] (A). 
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Fig. 3 (color online). Q^F^'"''" (Q^) as a function of . Left panel shows "benchmark theoretical data" 
(open circles) that include systematic uncertainties from different sources, explained in the text, in comparison 
with real experimental data Right panel serves to show the large form- factor behavior using a 

logarithmic plot. The shaded (green) strip contains the results for the BMS bunch of tt DAs, with the BMS 
model 18] being represented by a solid line. Also shown are the predictions of Agaev et al. [13] for their models 
I and III, extrapolated in terms of dashed (blue) lines to the remote timelike point 112 GeV^ from [211] . 

Figure [3] shows the calculated pion-photon transition form factor ^5] in comparison with various sets 
of data, specified in the caption. The graphics in the left panel displays the theoretical results in the 
form of "data" at the same momentum-transfer values as the experimental data and with theoretical 
uncertainties included as error bars — see Table I in the first work of Ref. [5j] for the numerics. The right 
panel serves to effect the large-Q^ behavior of the scaled form factor — therefore, a logarithmic scale. 
The shaded green strip embodies the predictions obtained in as described above. One observes: (i) 
The shaded strip fits the CELLO, CLEO and the 'CLEO regime' set of the BABAR data very precisely 
(overall < !)■ (ii) The prediction approaches asymptotically the limit \f2j.^ from below and agrees 
very well with the BABAR data for the 'y*'y'r]{r]') transition form factor, being in conflict with the steep 
rise of the BABAR data for Q^F'<''^'^° {Q"^) above 9 GeV^. (iii) The LCSR predictions of ^ (single 
solid blue lines) are not really reproducing the behavior of the BABAR data, despite the opposite claims 
by the authors and the use of pion DAs with a high number of harmonics. In fact, the reinforcement 
of the form factor can only be achieved by an enhancement of the endpoint region a; = 0, 1 of the 
pion DA, with the best agreement being provided by the flattop DA which includes all x values in an 
equal amount, inevitably entailing a worse fit to the data in the 'CLEO regime' and below (see [51 for 
details) . 

4 Conclusions 

To conclude, we cannot predict the rise of the scaled pion-photon transition form factor, observed by 
the BABAR Collaboration, using the standard scheme of QCD. Nor can we explain it in hindsight within 
the same context — despite opposite claims in jT3] and elsewhere [12]. Therefore, the BABAR data are 
unpredictable in QCD, but they are not unexplainable — in particular contexts which, however, are 
contingent on assumptions that can hardly be accommodated within the standard scheme of QCD. 
Moreover, there is a dichotomy of the pseudoscalar meson sector into two very disparate groups, one 
obeying the QCD asymptotics, the other not, depending on the particular shape of the pion DA. Indeed, 
for BMS-like models (double- humped and endpoint-suppressed) the predictions completely agree 
with the CELLO [H and all data in the 'CLEO region ' nil for the pion-photon transition form factor, 
as well as with the BABAR data 6] for the 77 ( 77')-photon transition form factor, while being in strong 
disagreement with the BABAR data ^] for the Tr'^-photon form factor above 9 GeV^. On the other hand, 
using a flatlike DA [H; [13; [1^ (which describes an unrealistic pointlike pion) reproduces the 7*77r" 
BABAR data but fails to comply with those on the 77 ( 77')-photon transition leading to predictions above 
the experimental data [5|] (see also [1^ [13; [13 )• These findings indicate a possible strong violation 
of the isospin symmetry in the pseudoscalar meson sector {tt° and ris) that has not been observed in 
other measurements so far. Therefore, the reproducibility of the BABAR results for the pion-photon 
transition form factor by other experiments (e.g., by Belle) is of paramount importance. 
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